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Abstract. In this paper, I present and discuss some recent observational 
results which may have important implications for our understanding of late 
phases of binary evolution. 
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1. INTRODUCTION 

Symbiotic stars are interacting binaries, consisting of an evolved giant (either a 
normal red giant in S-types or a Mira variable embedded in an optically thick dust 
shell in D-type) transferring mass to a hot and luminous white dwarf (although in 
a few cases a neutron star has been found). Since the binary must have enough 
room to accommodate the red giant, and in the case of D- types also its dust shell, 
symbiotic stars have the longest orbital periods and so the largest component 
separations among interacting binaries. The interacting stars are surrouned by a 
rich and luminous circumstellar environment resulting from the presence of both 
an evolved giant with a heavy mass loss and of a hot companion abundant in 
ionizing photons and often emanating its own wind. In particular, strongly differ- 
ent surroundings are expected, including both ionized and neutral regions, dust 
forming regions, accretion and excretion disks, interacting winds and jets. Such a 
complex structure makes symbiotic stars very important tracers of late stages of 
stellar evolution of low and intermediate-mass stars, and an excellent laboratory 
for exploring various aspects of interactions and evolution in binary systems (e.g. 
Corradi, Mikolajewska & Mahoney 2003, and references therein). Firming links 
between symbiotic stars and other families of interacting binaries with evolved 
giants is essential to understand, for instance, the role of such binaries in the for- 
mation of stellar jets, planetary nebulae, novae, super soft X-ray sources (SSXS), 
and Type la Supernovae (SNe la). Many of these are issues concerning the late 
stages of stellar evolution which are presently poorly known, but with important 
implications on our understanding of the stellar population and chemical evolution 
of galaxies, as well as extragalactic distance scale. 

Although significant progress has been done in deriving physical parameters 
of symbiotic stars, including orbital parameters, radii and mass loss rates (see 
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Fig. 1. The orbital period distribution. The shaded regions denote different popula- 
tions (upper panel) - the MCs and yellow, disk and buldge, as well as different activity 
classes (lower panel) - recurrent and ordinary symbiotic novae, systems with Z And-like 
activity, and systems without recorded outbursts. 
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Mikolajewska 2007, 2008, 2011a for reviews), there are important issues concerning 
their interactions and evolution which are not yet fully understood. 

The aim of this paper is to address some problems posed by observations of 
symbiotic binaries to theoretical models and binary evolution scenarios. 

2. ORBITAL PERIOD DISTRIBUTION 

Significant progress has been also made in deriving orbital periods using mas- 
sive photometric survey databases (Gromadzki et al. 2011; Mikolajewska 2011a, 
and references therein) . They are now known for more than 50 % of all S-type sys- 
tems included in Belczyhski et al.'s catalog (2000), although many of these periods 
have not yet been confirmed by radial velocity studies. Gromadzki et al. (2011) 
have argued that this makes ~ 80% of orbital periods that can be determined 
from photometric observations which means that the resulting period distribution 
should not be affected by selection effects. 

Significant progress has been also made in deriving orbital periods using mas- 
sive photometric survey databases (Gromadzki et al. 2011; Mikolajewska 2011a, 
and references therein) . They are now known for more than 50 % of all S-type sys- 
tems included in Belczyhski et al.'s catalog (2000), although many of these periods 
have not yet been confirmed by radial velocity studies. Gromadzki et al. (2011) 
have argued that this makes ~ 80% of orbital periods that can be determined 
from photometric observations which means that the resulting period distribution 
should not be affected by selection effects. 

The distribution of orbital periods is shown in Fig. 1. First of all, despite 
of continuously increasing number of measured periods, the main characteris- 
tics of their distribution remain practically the same as in earlier studies (e.g. 
Mikolajewska 2003, 2007). In particular, they peak around 600 d, and less than 
30 % systems have P or b above 1000 d. One of extremely important issue addressed 
by these results is that the currently used population synthesis models (PSM) fail 
to reproduce the observed orbital period distribution. The PSM predicts a P or b 
distribution peaked at ~ 1500 d and only 20% symbiotic systems with periods 
below 1000 d (Lu et al. 2006). The key to solve this problem is more advanced 
approach to mass transfer in symbiotic stars, and in fact in any interacting binaries 
involving red giants (e.g. Podsiadlowski & Mohamed 2007). In fact, there is clear 
evidence that in majority of S-type systems mass transfer and strong tidal interac- 
tion must have taken place before the present white dwarf was formed. Practically 
all systems with P or b below ~ 900 d have circular orbits, and for those with white 
dwarf mass above ~ 0.5 M , the present binary separation is too small and the 
orbital period too short to accommodate the former AGB star. 

Fig. 1 shows that there is no significant difference between the orbital period 
distribution for the symbiotic systems in the Galactic disk and those in the buldge. 
The yellow systems with a G-K giant donor and those in the MCs may split into 
two populations, with P or b ~ 200-600 d (only the low metallicity halo systems 
with K-type giants) and > 1000 d (dominated by the MC systems and D'-type 
systems), respectively. 

Fig. 1 also shows the P or b distribution for different activity classes. It is 
remarkable that vast majority (~ 90 %) of symbiotic stars with multiple outburst 
Z And-type activity have orbital periods below ~ 1000 d, whereas all four symbiotic 
recurrent novae (SyRNe) have orbital periods below 600 d. This may be related 
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Fig. 2. Left: Orbital periods as a function of the spectral types. Different symbols 
represent different activity classes: recurrent (stars) and ordinary (diamonds) symbiotic 
novae, systems with Z And-like activity (triangles), systems without recorded outbursts 
(circles) and systems with enhanced activity related to the periastron passage (squares) . 
The filled symbols denote objects included in Miirset & Schmid (1999). The lines make 
the period limits found by Miirset & Schmid (1999). Right: The Roche lobe radius versus 
the cool giant's radius, as derived from the spectral type, for systems with K and M-type 
giants. Lines denote -Rrl = R, -Rrl = 2R, and _Rrl = 4 7?, respectively. 

to higher mass transfer and accretion rates in the SyRNe and Z And-like systems 
than in remaining systems. 

3. MASS TRANSFER AND ACTIVITY 

One of the most fundamental question concerning S-type symbiotic stars is the 
mode of mass transfer, via stellar wind or Roche lobe overfow (RLOF), and the 
possibility of an accretion disk formation. 

Among the main arguments against the RLOF are the red giant radii derived 
from their rotational velocities, which, in most cases indicate R g ~ 0.5 i?RL, and 
the relation for the minimum orbital period for symbiotic systems with red giants 
of a given spectral type found by Miirset & Schmid (1999) which they interpreted 
as a restriction to the radius of the red giant. They also argued that this limit 
corresponds to a configuration where the radius of the giant reaches about half 
distance from its center to the inner Lagrangian point, L\. However, their relation 
seems to be not so obvious with the present sample of systems with known periods, 
a factor of 3 larger than that available for Miirset & Schmid. The orbital periods 
versus the spectral types together with the period limits found by Miirset & Schmid 
are plotted in the left panel of Fig. 2, whereas the right panel shows the Roche lobe 
radius (which, in our opinion, is more relevant that the distance to L\) versus the 
average radius for the corresponding spectral type for systems with K and M-type 
giants. The average Roche lobe radius (Paczyhski 1971) was calculated adopting 
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the mass ratio, M g /M w( j = 2.5, and the total mass of 2M (Mikolajewska 2003). 
Many systems appear to have orbital periods significantly shorter than the lower 
limit of Miirset & Schmid, and their spectral types indicate R g larger (and even 

close to) i?RL- 

Moreover, there is clear evi- 
dence for ellipsoidal light curve 
variations in many S-type systems 
with P or b below ~ 1000 d, in- 
cluding three of the four known 
SyRNe (Mikolajewska 2003, 2007, 
2011a; Rutkowski et al. 2007; 
Schaefer 2009; Yudin et al. 2005). 
Such changes can be, however, de- 
tected only in systems with rel- 
atively high inclination, and, in 
addition, they usually show up 
only in red and near infrared light 
(where the red giant dominates), 
whereas for most symbiotic stars 
only visual light curves are avail- 
able (Fig. 3). This is particu- 
larly unfortunate in the case of ac- 
tive Z And-like systems and inac- 
tive systems with very luminous 
hot component like SY Mus which 
dominate among S-type symbi- 
otics with shorter periods. Any- 
way, it seems that ellipsoidal vari- 
ability is present in IJHK light 
curves of all multiple outburst Z 
And-type systems, and in some 
steady systems as well, with the 
orbital inclination above <~ 60° 
(Rutkowski et al. 2007). Such 
variability is, however, definitely 
absent in the symbiotic nova AG 
Peg (Rutkowski et al. 2007). 

The presence of tidally distorted giants in symbiotic systems would facilitate 
accretion disk formation. This is particularly interesting in relation to the observed 
outburst characteristics of Z And-like systems which can be qualitatively accounted 
for by unstable accretion disk around the thermonuclear H-shell burning white 
dwarf (Mikolajewska 2003; Sokoloski et al. 2006) and the symbiotic recurrent 
novae where high accretion rates required to account for their short recurrence 
timescales could not be easily achieved by wind accretion. There is also continual 
problem with interpretation of spectra of some symbiotics, including RW Hya and 
SY Mus, the stable systems with tidally distorted giants. Their far ultraviolet 
spectral energy distributions indicate a factor of ~ 2 lower temperatures than 
emission lines which may indicate the presence of accretion disks also in these 
systems (e.g. Sion et al. 2002; Kolb et al. 2004; Sion et al. 2008). 
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Fig. 3. V (ASAS database; Pojmariski 2002) 
and J (Rutkowski et al. 2007) light curves of SY 
Mus phased with the orbital period of 624.5 d. 
The ellipsoidal variability is visible only in the 
near infrared light. 
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There are, however, two embarrassing problems with the ellipsoidal variability 
of the symbiotic giants. 

First, the mass ratio, measured by radial velocity curves for systems like AR 
Pav, CI Cyg and BF Cyg, is of ~ 2-4 (where the red giant mass donor is the more 
massive component), and according to the binary evolution theory they should 
undergo dynamically unstable mass transfer. One may argue that all these systems 
are about to enter unstable mass transfer, however, the very short timescales for 
this phase (a few 10 3 yrs or so) make it very unlikely. Thus far efforts to find new 
mechanisms to stabilize mass transfer have not been successful (e.g. Podsiadlowski 
& Mohamed 2007; Chan et al. 2010). 

Second problem faced by most symbiotic systems with ellipsoidal variability 
is the systematic discrepancy between the radii derived from light curve analysis 
and those measured from rotational velocities (Mikolajewska 2007). The only 
systems which do not show such an effect are CI Cyg and AX Per. To make the 
situation even more complicated, most researchers believe that the giant donor in 
the symbiotic recurrent nova RS Oph underfills its Roche lobe by a factor of <~ 2 
whereas its rotational velocity is consistent with a synchronized Roche lobe-filling 
star (Brandi et al. 2009). 

There is no obvious solution to the 'radius-u rot sini problem in tidally inter- 
acting binaries. Asynchronous rotation seems unlikely because of circular orbits 
and short timescales for the synchronization (actually shorter than that for cir- 
cularization) in such binaries. Although it is possible that the assumption of a 
spherical giant with a simple limb-darkening biases the v ro t sini measurements 
(e.g. Orosz & Hauschildt 2000), the observed discrepancy seems to be too large 
for such an explanation. It is also hard to believe that the ellipsoidal variations 
are caused by the very extended stratified atmospheres (with molecular layers) or 
slow stellar wind filling the Roche lobe rather than the star itself because the effect 
is independent of the spectral type (from M2 on RW Hya to M6 in AR Pav), and 
the multicolor light curves can be fitted with single atmosphere ellipsoidal model 
(e.g. Yudin et al. 2005; Rutkowski et al. 2007). At present the most promising 
explanation seems to be modification of the Roche potential used for the ellipsoidal 
variability by including radiation effects and pulsations. For example, Schuerman 
(1972) analytically showed that the contact surface can shrink up to 2-3 times 
for stars with strong winds. Later extensions of these calculations showed that 
symbiotic stars can be good candidates for such an adapted Roche potential (Der- 
mine et al. 2009). In fact, the cool giants in most S-type systems (including SY 
Mus) show low amplitude pulsations (Gromadzki et al. 2007, 2011) whereas the 
radio observations indicate strong, ~ 10~ 7 M Q /yr, winds (e.g. Mikolajewska et al. 
2002). 

All these possibilities can be tested by direct measurements of the radii of the 
red giants in brightest and closest systems by interferometric technique. After all, 
it may happen that synchronization model for binaries also needs improvement, 
for instance, by including stellar wind and pulsations. 

4. SYMBIOTIC STARS AND PROGENITORS OF TYPE IA SUPERNOVAE 

Understanding symbiotic stars may also help to solve one of the timeliest prob- 
lems in modern astrophysics - missing progenitors of SNe la. There is general con- 
sensus that they result from thermonuclear disruption of CO white dwarf reaching 
the Chandrasekhar mass either due to mass accretion from a non-degenerated com- 
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panion (single-degenerate or SD model) or mass transfer between and/or merger 
of two white dwarfs (double-degenerate or DD model). Unfortunately, the pro- 
genitors of SNe la have never been observed before the explosion, and each of the 
proposed scenarios has its pros and cons (DiStefano & Orio 2011, and references 
therein). With regard to this problem symbiotic stars are important for several 
reasons. First, recent observations strongly support the single-degenerate model 
for at least a subset of SNe la, in particular systems with red-giant companions re- 
sembling the S-type symbiotic RS Oph (Patat et al. 2011, and references therein). 
Second, as recently pointed out by DiStefano (2010), the key to solve the missing 
progenitor problem lies in understanding the appearance of nuclear-burning white 
dwarfs which according to theoretical predictions should appear as bright SSXSs 
whereas observations show that they are too few if them in all types of galax- 
ies to serve as progenitors of SNe la, and moreover it applies to in any scenario. 
Whether or not symbiotic stars can be SNe la progenitors, most of them contain 
nuclear-burning white dwarfs (e.g. Mikolajewska 2003, 2011a) but the only SSXSs 
associated with symbiotic stars are two systems in the SMC, SMC 3 and LIN 358, 
and AG Dra, the galactic system with the lowest measured metallicity. Moreover, 
symbiotic stars appear to be promising candidates for both SD and DD scenarios 
because some systems like AR Pav, if they will pass through a common envelope 
phase, may become close double white dwarf systems with total mass comparable 
to the Chandrasekhar limit (Mikolajewska 2011b). 

5. CONCLUSIONS 

Summarizing, symbiotic binaries continually pose challenges to the binary evo- 
lution theory. Understanding the mass transfer and accretion processes in these 
systems is not only essential for understanding their present status but also for 
estimating their rates as well as for any interacting binaries involving evolved 
giants. 
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